REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  NO.  0704-0188 


The  public  reporting  burden  for  this  coilection  of  information  is  estimated  to  average  1  hour  per  response,  inciuding  the  time  for  reviewing  instructions, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  compieting  and  reviewing  the  coiiection  of  information.  Send  comments 
regarding  this  burden  estimate  or  any  other  aspect  of  this  coilection  of  information,  including  suggesstions  for  reducing  this  burden,  to  Washington 
Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Ariington  VA,  22202-4302. 
Respondents  shouid  be  aware  that  notwithstanding  any  other  provision  of  iaw,  no  person  shaii  be  subject  to  any  oenaity  for  failing  to  comply  with  a  coiiection 
of  information  if  it  does  not  dispiay  a  currentiy  vaiid  OMB  controi  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


2.  REPORT  TYPE 
Final  Report 


1.  REPORT  DATE  (DD-MM-YYYY) 
11-02-2014 


4.  TITLE  AND  SUBTITLE 

Unsteady  Aerodynamics  of  Flapping  Wings  at  Re=10.000  - 
100,000  for  Micro-Air  Vehicles-Final  Report 


3.  DATES  COVERED  (From  -  To) 
1-Nov-2008  -  31-Oct-2013 


5a.  CONTRACT  NUMBER 
W911NF-08- 1-0433 


5b.  GRANT  NUMBER 


6.  AUTHORS 

Danesh  Tafti,  Kamal  Viswanath,  Venkata  Ravi  Kasibhotia 


5c.  PROGRAM  ELEMENT  NUMBER 
611102 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAMES  AND  ADDRESSES 

Virginia  Polytechnic  Institute  &  State  Unive 
North  End  Center,  Suite  4200 
300  Turner  Street,  NW 

Blacksburg,  VA  24061  -0001 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS 
(ES) 

U.S.  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 
ARO 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

53417-EG.12 


12.  DISTRIBUTION  AVAILIBILITY  STATEMENT 
Approved  for  Public  Release;  Distribution  Unlimited 


13.  SUPPLEMENTARY  NOTES 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  contmed  as  an  official  Department 
of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


14.  ABSTRACT 

The  overall  project  goal  was  to  investigate  low  Reynolds  number  (10,000-100,000)  flapping  flight  for  MAV 
applications.  Large-Eddy  Simulations  (EES)  is  used  for  unsteady  aerodynamic  analysis.  The  dominant  force 
producing  mechanism  is  the  formation  of  a  spiral  leading-edge-vortex  (LEV).  Effect  of  wind  gusts  and  stroke 
deviation  is  investigated  followed  by  the  aerodynamics  of  a  flapping  fruit  bat  wing.  Einally,  very  high  Reynolds 
number  dynamic  stall  on  a  pitching  NACA-012  airfoil  is  studied.  The  work  resulted  in  3  archived  journal 

_ — - 1  nu  rv  — a  1  a/t  c  „aau:^ i  - 1 _ — „  — „ — a  — 


15.  SUBJECT  TERMS 

Flapping  flight.  Micro-air  vehicles.  Bat  flight.  Dynamic  stall 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

15.  NUMBER 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF  PAGES 

UU 

UU 

UU 

UU 

Danesh  Tafti 


19b.  TELEPHONE  NUMBER 
540-231-9975 


Standard  Form  298  (Rev  8/98) 
Prescribed  by  ANSI  Std.  Z39. 18 


Report  Title 

Unsteady  Aerodynamics  of  Flapping  Wings  at  Re=10.000  -  100,000  for  Micro-Air  Vehicles-Final  Report 

ABSTRACT 

The  overall  project  goal  was  to  investigate  low  Reynolds  number  (10,000-100,000)  flapping  flight  for  MAV  applications.  Large-Eddy 
Simulations  (LES)  is  used  for  unsteady  aerodynamic  analysis.  The  dominant  force  producing  mechanism  is  the  formation  of  a  spiral  leading- 
edge-vortex  (LEV).  Effect  of  wind  gusts  and  stroke  deviation  is  investigated  followed  by  the  aerodynamics  of  a  flapping  fruit  bat  wing. 
Finally,  very  high  Reynolds  number  dynamic  stall  on  a  pitching  NACA-012  airfoil  is  studied.  The  work  resulted  in  3  archived  journal 
publications,numerous  conference  proceedings,!  Ph.D.  and  1  M.S.  thesis.  Two  additional  journal  papers  are  under  preparation. 


Enter  List  of  papers  submitted  or  published  that  acknowledge  ARO  support  from  the  start  of 
the  project  to  the  date  of  this  printing.  List  the  papers,  including  journal  references,  in  the 
following  categories: 

(a)  Papers  published  in  peer-reviewed  journals  (N/A  for  none) 


Received  Paper 


08/23/2013  7  00  Kamal  Viswanath,  Tafti  D.  K..  Effect  of  Stroke  Deviation  on  Forward  Flapping  Flight, 
AIAA  Journal,  (01  2013):  145.  doi: 

08/29/2011  1 .00  Kamal  Viswanath,  Danesh  K.  Tafti.  Effect  of  Frontal  Gusts  on  Forward  Flapping  Flight, 
AIAA  Journal,  (9  2010):  0.  doi:  10.2514/1.J050263 

TOTAL:  2 


Number  of  Papers  published  iu  peer-reviewed  jouruals: 


(b)  Papers  published  in  non-peer-reviewed  journals  (N/A  for  none) 


Received  Paper 


TOTAL: 


Number  of  Papers  published  iu  uou  peer-reviewed  jouruals: 


(c)  Presentations 

Tafti,  D.  K.,  Viswanath,  K.  and  K.  Nagendra,  Deconstructing  the  Essential  Elements  of  Bat  Flight,  APS  66th  Annual  Meeting,  Division  of 
Fluid  Dynamics,  24-26  November  2013,  Pittsburgh,  PA. 


Number  of  Presentations:  1.00 

Non  Peer-Reviewed  Conference  Proceeding  publications  (other  than  abstracts): 


Received  Paper 


08/23/2013  8.00  K.  Viswanath,  K.  Nagendra,  D.  K.  Tafti.  Climbing  Flight  of  a  Fruit  Bat  Deconstructed, 

AIAA  Science  and  Technology  Forum  and  Exposition  2014.  13-JAN-14,  .  :  , 

08/29/2011  2.00  Danesh  Tafti.  UNSTEADY  AERODYNAMICS  OF  A  FLAPPING  WING  FOR  MICO  AIR  VEHICLE  (MAV) 
APPLICATIONS  , 

Proceedings  of  the  37th  National  &  4th  International  Conference  on  Fluid  Mechanics  and  Fluid  Power,  NT 
Madras,  Chennai,  India.  16-DEC-10,  .  :  , 

TOTAL:  2 


Number  of  Non  Peer-Reviewed  Conference  Proceeding  publications  (other  than  abstracts): 

Peer-Reviewed  Conference  Proceeding  publications  (other  than  abstracts): 


Received  Paper 


02/11/2014  11.00  K.  Viswanath,  K.  Nagendra,  D.  Tafti.  Climbing  Flight  of  a  Fruit  Bat  Deconstructed, 

AIAA  SCITECH  2014.  13-JAN-14,  .  :  , 

08/29/201 1  3.00  Kamal  Viswathan,  Danesh  Tafti.  Effect  of  Stroke  Deviation  on  Forward  Flapping  Flight , 
50th  AIAA  Aerospace  Sciences  Meeting.  09-JAN-12,  .  :  , 


TOTAL: 


2 


Number  of  Peer-Reviewed  Conference  Proceeding  publications  (other  than  abstracts): 


(d)  Manuscripts 


Received  Paper 


08/23/2013  9.00  K.  Viswanath,  K.  Nagendra,  J.  Cotter,  M.  Frauenthal,  D.K.  Tafti.  Straight-Line  Climbing  Flight 
Aerodynamics  of  a  Fruit  Bat, 

Physics  of  Fluids  (05  2013) 

08/27/2012  4  00  Kamal  Viswanath,  Danesh  K.  Tafti.  Effect  of  Stroke  Deviation  on  Forward  Flapping  Flight, 
AIAA  Journal  (10  2011) 

TOTAL:  2 


Number  of  Manuscripts: 


Books 


Received  Paper 


TOTAL: 


Patents  Submitted 


Patents  Awarded 


Awards 

D.  Tafti  awarded  William  S.  Cross  Professorship,  2009. 

D.  Tafti  awarded  College  of  Engineering  Deans  Researeh  Award,  2012 


Graduate  Students 


NAME 

PERCENT  SUPPORTED  Discipline 

Kamal  Viswanath 

0.75 

Venkata  Ravi  Kasibhotia 

0.08 

FTE  Equivalent: 

0.83 

Total  Number: 

2 

Names  of  Post  Doctorates 


NAME 

PERCENT  SUPPORTED 

FTE  Equivalent: 

Total  Number: 

Names  of  Faculty  Supported 


NAME 

PERCENT  SUPPORTED  National  Academy  Member 

Danesh  Tafti 

0.06 

FTE  Equivalent: 

0.06 

Total  Number: 

1 

Names  of  Under  Graduate  students  supported 


NAME 

PERCENT  SUPPORTED  Discipline 

Jeffrey  Cotter 

0.00 

Mathew  Fruenthal 

0.00 

FTE  Equivalent: 

0.00 

Total  Number: 

2 

Student  Metrics 

This  section  only  applies  to  graduating  undergraduates  supported  by  this  agreement  in  this  reporting  period 

The  number  of  undergraduates  funded  by  this  agreement  who  graduated  during  this  period: .  0.00 

The  number  of  undergraduates  funded  by  this  agreement  who  graduated  during  this  period  with  a  degree  in 

seienee,  mathematics,  engineering,  or  technology  fields: . 0.00 

The  number  of  undergraduates  funded  by  your  agreement  who  graduated  during  this  period  and  will  continue 

to  pursue  a  graduate  or  Ph.D.  degree  in  science,  mathematics,  engineering,  or  technology  fields: . 0.00 

Number  of  graduating  undergraduates  who  achieved  a  3.5  GPA  to  4.0  (4.0  max  scale): . 0.00 

Number  of  graduating  undergraduates  funded  by  a  DoD  funded  Center  of  Excellence  grant  for 

Education,  Research  and  Engineering; . o.OO 

The  number  of  undergraduates  funded  by  your  agreement  who  graduated  during  this  period  and  intend  to  work 

for  the  Department  of  Defense . 0.00 

The  number  of  undergraduates  funded  by  your  agreement  who  graduated  during  this  period  and  will  receive 

scholarships  or  fellowships  for  further  studies  in  science,  mathematics,  engineering  or  technology  fields: . 0.00 


Names  of  Personnel  receiving  masters  degrees 

NAME 

Venkata  Ravi  Kasibhotfa 

Total  Number:  1 


Names  of  personnel  receiving  PHDs 


NAME 

Kamal  Viswanath 

Total  Number: 

1 

Names  of  other  research  staff 


NAME  PERCENT  SUPPORTED 

FTE  Equivalent: 

Total  Number: 

Sub  Contractors  (DD882) 

Inventions  (DD882) 


Scientific  Progress 


Considerable  scientific  progress  was  made  during  the  course  of  this  project.  Our  initial  focus  was  on  flat  plate  (stiff  and  flexible) 
flapping  wing  kinematics.  In  this  geometry  we  investigated  the  effect  of  large  gusts  engulfing  the  wing  and  the  effect  of  stroke 
deviations  on  the  aerodynamic  forces.  Then  we  used  the  knowledge  gathered  in  this  domain  to  attack  the  complex  measured 
kinematics  of  a  bat  wing.  The  experiments  were  done  at  Brown  University.  We  used  the  point  traces  on  the  wing  to  construct  a 
space-time  mesh  of  the  wing  surface  which  was  then  immersed  in  a  background  grid.  One  surprising  conclusion  from  this  work 
was  that  in  spite  of  the  apparent  complexity  of  wing  motion,  the  motion  could  be  decomposed  into  2-3  modes  which  reproduced 
the  force  signature  of  the  wing.  Finally,  we  investigated  the  challenging  problem  of  dynamic  stall  at  very  high  Reynolds  numbers 
on  a  pitching  NACA0012  airfoil.  We  found  that  the  upstroke  behaved  quite  benignly  but  the  downstroke  produced  very  complex 
flow  structures  and  the  lift  was  very  dependent  on  the  shedding  and  formation  of  LEVs  and  TEVs.  Fair  comparisons  were 
obtained  with  experiments  considering  that  the  experimental  tunnel  cross-section  is  much  smaller  than  the  computational 
domain  and  could  potentially  interfere  with  the  vorticity  dynamics  on  the  airfoil. 


Effect  of  Frontal  Gusts  on  Forward  Flapping  Flight: 

The  response  of  a  rigid  flapping  wing  in  forward  flight,  at  Re=1 0,000,  subjected  to  frontal  gusts  has  been  investigated.  The 
phasing  and  duration  of  the  gusts  and  their  impact  on  the  various  unsteady  mechanisms  are  analyzed  within  a  single  flapping 
cycle.  The  gust  is  characterized  by  a  step  function  with  integral  length  scale  much  larger  than  that  of  the  physical  dimension  of 
the  MAV  and  with  time  scale  much  smaller  than  the  flapping  time  period.  The  instantaneous  lift  and  thrust  profiles  were 
observed  to  be  influenced  by  a  combination  of  the  effective  angle  of  attack,  wing  rotation  and  the  Leading  Edge  Vortex  (LEV) 
structures  existing  in  the  flow  at  any  given  time,  with  the  LEVs  themselves  being  influenced  by  the  duration  and  magnitude  of 
the  change  in  effective  angle  of  attack.  Frontal  gusts  applied  during  the  downstroke  accelerated  the  development  of  the  flow 
resulting  in  the  formation  and  detachment  of  multiple  LEVs  on  the  wing  surface  which  increased  the  lift  and  thrust  illustrating  the 
importance  of  the  LEV  dynamics  to  force  production.  The  effect  of  the  gust  is  observed  to  be  diminished  when  it  occurs  during 
rapid  supination  of  the  wing.  The  lift  and  thrust  profiles  are  found  to  react  in  a  similar  fashion  for  gusts  applied  during  the 
downstroke  while  they  experienced  opposite  effects  during  the  upstroke.  During  the  upstroke,  force  characteristics  are  seen  to 
primarily  react  to  effective  AOA  changes  more  than  changes  in  flow  structures. 


Effect  of  Frontal  Gusts  on  Flexible  Wings  in  Forward  Flapping  Flight: 

The  response  of  a  flexible  wing  and  a  rigid  wing  in  flapping  forward  flight,  at  Re=1 0,000,  subjected  to  frontal  gusts  is 
investigated.  The  impact  of  a  gust  of  chosen  phasing  and  duration  on  the  various  unsteady  mechanisms  are  analyzed  within  a 
single  flapping  cycle.  The  gust  is  characterized  by  a  step  function  with  integral  length  scale  much  larger  than  that  of  the  physical 
dimension  of  the  MAV  and  with  time  scale  much  smaller  than  the  flapping  time  period.  The  instantaneous  lift  and  thrust  were 
observed  to  be  influenced  by  a  combination  of  the  effective  angle  of  attack,  chosen  kinematics,  transient  Leading  Edge  Vortex 
(LEV)  structures,  and  the  dynamic  cambering  in  the  case  of  the  membrane  wing.  The  effect  of  the  gust  on  the  lift  profile  is 
observed  to  have  opposite  effects  for  the  rigid  vs.  flexible  case  while  the  thrust  behaves  in  a  similar  profile. 

Effect  of  Stroke  Deviation  on  Forward  Flapping  Flight: 

The  performance  of  a  rigid  thin  surface  flat-plate  flapping  wing  in  forward  flight,  at  Re=1 0,000,  using  different  stroke  deviation 
trajectories  has  been  investigated  to  assess  the  different  capabilities  that  such  kinematics  might  offer.  The  instantaneous  lift  and 
thrust  profiles  were  observed  to  be  influenced  by  a  combination  of  the  Leading  Edge  Vortex  (LEV)  and  the  Trailing  Edge  Vortex 
(TEV)  structures  existing  in  the  flow  at  any  given  time.  Unlike  regular  no  deviation  flapping  cycles,  the  TEV  is  shown  to  be 
significant  for  out-of-plane  trajectories.  Both  clockwise  and  anticlockwise  variations  for  trajectory  choice  are  analyzed  for  their 
efficacy  as  an  improved  kinematic  choice  over  a  no-deviation  base  case.  The  power  requirements  for  the  different  cases,  based 
on  the  fluid  torques,  are  used  as  an  index  of  the  cost  of  performance  across  all  cases.  The  anti-clockwise  8-cycle  deviation  is 
shown  to  be  very  complex  with  high  power  costs  albeit  having  better  performance.  The  clockwise  0-cycle  holds  promise  in 
being  utilized  as  a  viable  stroke  deviation  trajectory  for  forward  flight  over  the  base  case. 

Straight-Line  Climbing  Flight  Aerodynamics  of  a  Fruit  Bat: 

From  flight  data  obtained  on  a  fruit  bat,  Cynopterus  brachyotis,  a  kinematic  model  for  straight-line  flapping  motion  is  extracted 
and  analyzed  in  a  CFD  framework  to  gain  insight  into  the  complexity  of  bat  <j,ight.  The  intricate  functional  mechanics  and 
architecture  of  the  bat  wings  set  it  apart  from  other  vertebrate  flight.  The  extracted  kinematic  model  is  simulated  for  a  range  of 
Reynolds  numbers,  to  observe  the  effect  these  phenomena  have  on  the  unsteady  transient  mechanisms  of  the  flow  produced 
by  the  flapping  wings.  The  Strouhal  number  calculated  from  the  data  is  high  indicating  that  the  oscillatory  motion  dominates  the 
flow  physics.  From  the  obtained  data,  the  bat  exhibits  fine  control  of  its  mechanics  by  actively  varying  wing  camber,  wing  area, 
torsional  rotation  of  the  wing,  forward  and  backward  translational  sweep  of  the  wing,  and  wing  conformation  to  dictate  the  fluid 
dynamics.  As  is  common  in  flapping  flight,  the  primary  force  generation  is  through  the  attached  unsteady  vortices  on  the  wing 
surface.  This  force  output  is  modulated  by  the  bat  through  varying  the  wing  camber  and  the  wing  area.  The  power  requirement 


for  the  kinematics  is  anaiyzed  and  correiated  with  the  aerodynamic  performance.  Understanding  the  physics  of  these  robust 
fliers  wiii  heip  devise  prescribed  kinematics  for  mechanicai  flappers  as  weii  as  improve  upon  them  from  nature. 


3D  Dynamic  Staii  Simuiation  of  NACA0012  Airfoii: 

The  work  investigates  the  physics  behind  the  dynamic  staii  process  by  simuiating  the  flow  past  pitching  NACA-0012  airfoii  at 
100,000  and  1,000,000  Reynoids  number  based  on  the  chord  iength  of  the  airfoii  and  at  different  reduced  frequencies  of  0.188 
and  0.25,  respectiveiy  in  a  three  dimensionai  flow  fieid.  The  mean  angies  of  attack  are  12  and  15  degrees  and  the  ampiitudes  of 
pitching  are  6  and  10  degrees,  respectiveiy.  The  turbuience  in  the  flow  fieid  is  resoived  using  iarge  eddy  simuiations  with 
dynamic  Smagorinsky  modei  at  the  sub  grid  scaie.  The  iift  hysteresis  piots  of  this  simuiation  for  both  the  configurations  are 
compared  with  the  corresponding  experiments.  The  deveiopment  of  dynamic  staii  vortex,  vortex  shedding  and  reattachment  as 
predicted  by  the  present  study  are  discussed  in  detaii.  There  is  a  fairiy  good  match  between  the  predicted  and  experimentaiiy 
measured  iift  coefficient  during  the  upstroke  for  both  cases.  The  net  iift  coefficient  for  the  Re=1 00,000  case  during  downstroke 
matches  with  the  corresponding  experimentai  data,  the  present  study  under-predicts  the  iift  coefficient  as  compared  to  the 
experimentai  vaiues  at  the  start  of  downstroke  and  over-estimates  for  the  remaining  part  of  the  downstroke.  The  trend  of  the  iift 
coefficient  hysteresis  piot  with  the  experimentai  data  for 

the  Re=1 ,000,000  case  is  aiso  simiiar.  This  present  simuiations  have  shown  that  the  downstroke  phase  of  the  pitching  motion  is 
strongiy  three  dimensionai  and  is  highiy  compiex,  whereas  the  flow  is  practicaiiy  two  dimensionai  during  the  upstroke. 
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